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Abstract 

Possible modification in the velocity distribution in the non-resonant reaction rates leads 
to an extended reaction rate probability integral. The closed form representation for these 
thermonuclear functions are used to obtain the stellar luminosity and neutrino emission rates. 
The composite parameter C that determines the standard nuclear reaction rate through the 
Maxwell-Boltzmann energy distribution is extended to C* by the extended reaction rates 
through a more general distribution than the Maxwell-Boltzmann distribution. The new dis- 
tribution is obtained by the pathway model introduced by Mathai in 2005 [Linear Algebra 
and Its Applications, 396, 317-328]. Simple analytic models considered by various authors 
are utilized for evaluating stellar luminosity and neutrino emission rates and are obtained in 
generalized special functions such as Meijer's G-function and Fox's i?-function. The standard 
and extended non-resonant thermonuclear functions are compared by plotting them. Behav- 
ior of the new energy distribution, more general than Maxwell-Boltzmann is also studied. 

Keywords: Thermonuclear function, pathway model, reaction rate probability integral, stel- 
lar model, fusion energy, G-function. 

1 Introduction 

In one way or another, the mystery behind the distant universe is explored by the 
understanding of the Sun, the star nearest to Earth. It is the only star whose mass, 
radius and luminosity are known to highest accuracy. Structural change in the Sun is 
due to the evolution of the central gravitationally stabilized solar fusion reactor. Solar 
nuclear energy generation and solar neutrino emission are governed by chains of nu- 
clear reactions in this gravitationally stabilized solar fusion reactor [H [10] . Qualitative 
calculations of specific thermonuclear reaction rates require a large amount of experi- 
mental input and theoretical assumption. By using the theories from nuclear physics 
and kinetic theory of gases one can determine the reaction rate for low-energy non- 
resonant thermonuclear reactions in a non-degenerate plasma [9]. The formalization 
of the calculation of the reaction rate of interacting particles under cosmological or 
stellar conditions were presented by many authors [TJ [21]. For the most common case, 
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a nuclear reaction in which a particle of type 1 strikes a particle of type 2 producing a 
nucleus 3 and a new particle 4 is symbolically represented as 



l + 2^3 + 4 + £ 12 , (1.1) 

where is the energy release given by E%2 = (mi + rri2 — 7713 — m^c 2 , where rrii,i = 
1, 2, 3, 4 denote the masses of the particles, c denotes the velocity of light. The reaction 
rate r 12 of the interacting particles 1 and 2 is obtained by averaging the reaction cross 
section over the normalized density function of the relative velocity of the particles [T3"| 
[211 16]. Let ri\ and n 2 denote the number densities of the particles 1 and 2 respectively 
and cr(v) be the reaction cross section, where v is the relative velocity of the particles 
and f(v) is the normalized velocity density, then the reaction rate 7*12 is given by 

j-12 = (l - ^12 J n 1 n 2 (crf)i2 = (l - 7^12 J nin 2 ^ w(v)/(i;)dt> 

1 \ f°° /2E\ * 

1 - 2*12 J "ma 1 <r(E) f(E)dE, (1.2) 

where #12 is the Kronecker delta which is introduced to avoid double counting in the 
reaction, if particles 1 and 2 are identical. The quantity (crv)i2 is the thermally averaged 
product which is in fact the probability per unit time that two particles 1 and 2, confined 
to a unit volume, will react with each other. Here, /1 is the reduced mass of the particles 
given by \x = , E = is the kinetic energy of the particles in the center of mass 

system. From the literature [SI IZ3] , it may be noted that all the analytic expressions 
for astrophysically relevant nuclear reaction rates underline the hypothesis that the 
distribution of the relative velocities of the reacting particles always remains Maxwell- 
Boltzmannian for a non-relativistic non-degenerate plasma of nuclei in thermodynamic 
equilibrium. The Maxwell-Boltzmann relative kinetic energy distribution can be written 
as 

f M B D (E)dE = 2vr f-j-) 2 exp f-^j VEdE, (1.3) 
where k is the Boltzmann constant and T is the temperature. Substituting ( 1.3[ ) in 



(1.2) we get, 



1 \ f 8 Y / 1 \^ f 00 ^ / ( E\ 



r l2=[ l-- S ^ ni n.^-j j o MB )exp^-jd^. (1.4) 

The thermonuclear fusion depends on three physical variables, the temperature T, the 
Gamow energy Eg, and the nuclear fusion factor S(E). If two nuclei of charges Z\e 
and Z 2 e collide at low energies below the Coulomb barrier, the Gamow energy Eq is 
given by (26J, [1] 

E G = 2fx{7iaZ 1 Z 2 c) 2 , (1.5) 
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where a is the electromagnetic fine structure constant given by 



a 



e 

he 



1.6) 



where e is the quantum of electric charge, h is Planck's quantum of action, and a 
is ]^ Thus the Gamow factor, which is determined by the electromagnetic force, 
and the nuclear fusion factor S(E) determine the nuclear reaction cross section at low 
energies for non-resonant charged particles as [21 [6] 



Eg 
E 



;i-7) 



and S(E) is the cross section factor which is often found to be constant or a slowly 
varying function of energy over a limited range of energy given by [211 Ej 



S(E)~S(0) + — E+-^- 



E = — M£ 



Substituting (|1.7|) and (|1.8|) in (|1.4|) we obtain, 

1 



1 - ~5i2 n x n 2 

2 J \TljJi 



1 

kf 



v=0 



f 2 



S'M(O) 



y=0 



x / E^exp 
'o 



'fcT 



E 



dE. 



1.9) 



This is the non-resonant reaction rate probability integral in the Maxwell-Boltzmann 
case. The closed form evaluation of this integral can be seen in a series of papers 
by Mathai and Haubold, see for example Haubold and Mathai [121 H3], Mathai and 
Haubold [21] • The main aim of the present work is to extend the reaction rate prob- 



ability integral given in (1.9) by replacing the Maxwell-Boltzmann energy distribution 



by a more general energy distribution called the pathway energy distribution obtained 
by using the pathway model of Mathai, introduced in 2005. 



The paper is organized as follows: In the next section we discuss a more general 
energy distribution than Maxwell-Boltzmann and obtain the extended reaction rate 
probability integral in the non-resonant case. We take advantage of the closed form 
representation of the extended thermonuclear reaction rate for finding the luminosity 
and the neutrino emission rate of a so-called non-linear stellar model under consideration 
in section 3. Section 4 is devoted to find the desired connection between stellar structure 
parameters and the neutrino emission of the stellar model by using the closed form 
analytic representation of the extended reaction rate. A comparison of the Maxwell- 
Boltzmann energy distribution with the pathway energy distribution is done with the 
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help of graphs in section 5. Also we try to discriminate the standard and extended 
reaction rates. Concluding remarks are included in section 6. 

2 Extended non-resonant thermonuclear reaction rate 
and its closed forms 

In recent years, possible deviations of the velocity distribution of the plasma particles 
from the Maxwell-Boltzmann in connection with the production of neutrinos in the 
gravitationally stabilized solar fusion reactor has been pointed out [21 El ESI EZ1 ESI ED] • 
It was initiated by Tsallis, the originator of non-extensive statistical mechanics [281 EH1 
[8], who has used the g-exponential function as the fundamental distribution instead 
of the Maxwell-Boltzmann distribution. An initial attempt to extend the standard 
theories of reaction rates to Tsallis statistics was done by many authors, see Mathai 
and Haubold [22], Saxena et al. [21] • In 2005, Mathai introduced the pathway model by 
which even more general distributions can be incorporated in the theory of reaction rates 
[2Ql [22] . Initially, the pathway model was introduced for the matrix variate case to cover 
many of the matrix variate statistical densities. The scalar case is a particular one there. 
Subsequently, Mathai, his co-workers and others found connection of the pathway model 
with information theory, fractional calculus, and Mittag-Leffler functions. The pathway 
model can be effectively used in any situation in which we need to switch between three 
different functional forms, namely generalized type-1 beta form, generalized type-2 
beta form, and generalized gamma form, using the pathway parameter q. For practical 
purposes of fitting experimental data, the pathway model can be utilized to switch 
between different parametric families with thicker or thinner tail. The pathway model 
for the real scalar case can be summarized as follows: 

f x [ x ) = cis 7_1 [l- a{l-q)x s )^, a>0,5>0,l-a(l-g)s 5 >0,7>0,g< 1 (2.1) 

is the generalized type-1 beta form of the pathway model. This is a model with right 
tail cut-off for q < 1. The Tsallis statistics for q < 1 can be obtained from this model 
by putting a = 1, 7 = 1, 6 — 1 [El EHJ [29]. Other cases available are the regular type-1 
beta density, Pareto density, power function, triangular and related models [23]. The 
generalized type-2 beta-form of the pathway model is given by 

f 2 (x) =c 2 x 7 - 1 [l + a(q-l)x 5 )~^, < x < 00, q > 1, 7 > 0, 8 > 0. (2.2) 

Here, also for 7 = 1, a = 1, 5 = 1, we get the Tsallis statistics for q > 1 [El EEJ 129] . 
Other standard distributions coming from this model are regular type-2 beta density, 
F-distribution, Levy model and related models [23]. When q — > 1, fi(x) and f2(x) will 
reduce to the generalized gamma form of the pathway model,,,,, given by 

/ 3 (x) = c 3 ^- 1 e-^,x>0. (2.3) 



This model covers generalized gamma, gamma, exponential, chisquare, Weibull, Maxwell- 
Boltzmann, Rayleigh and related densities. Here, ci, c 2 and c 3 are defined in (2.1 ), (2.2) 



and (2.3), respectively, and are the normalizing constants if we consider statistical den- 
sities. 

By a suitable modification of the Maxwell-Boltzmann distribution, given in (1.3), 
through the pathway model, we get a more general energy distribution called the path- 
way energy distribution given by the density 



f PD (E)dE 



27r(g- 1)1 r (i=i 



{nkT) 
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l + (g-r 



E 
kf 



i 

9-1 



dE, 



(2.4) 



for q > 1, — | > 0. The Maxwell-Boltzmann energy distribution can be retrieved 
from (2.4) by taking q — > 1. Thus the reaction rate probability integral given in (1.4) 
can be modified by using (2.4) and we get the extended reaction rate as 

r 

7T/i j \ kT 
E 



ri2 



1 - 2^12 ) n l n 2 



1 

9-1 



E 



S'M(O) 



u=0 



E v 



i + (g 



i 

"9-1 



exp 



Eg 
E 



dE 



(2.5) 



for q > 1, -^j- — I > 0. Substituting y = ^ and 2 = (fifr) 2 we obtain the above integral 



in a more convenient form 

1 



1 



812 n x n 2 



(9-1) 



3 V g-l 



E 

i/=0 



1 \ 2 
kT 



S'W(O) 



+ (g - l)y] ^e" 



for g > 1, ^-j- — I > 0. Here we consider the integral to be evaluated as 



y v [l + (g - l)y] 



-xy 



'dy. 



(2.6) 



(2.7) 



The integral can be evaluated by the techniques in applied analysis and can be obtained 
in closed form via Meijer's G-function as PUJ [TTJ US] 



7T 2 



llq 



^1, 



,1 / (g- i)^ 2 |2-^r- 



4 



I0,i,f+1 



(g - i)^+ir ^ "'" 1 ' (2 ' 8) 

which yields the non-resonant reaction rate probability integral in the extended 



1 - -^12 ) ™1™2 



7T 



-1 



E 

i/=0 



g — 1 N 
~kf 



x 



S H (0)„ 3 ,1 
^1,3 



(q-l)E G ,2--^- 



AkT 



I o,|^+i 



(2.9) 
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Meijer's G-function and its properties can be seen in Mathai and Saxena [21], Mathai 
[19J. We can obtain series expansions of the G-function given in (2.9) by combining the 
theories of residue calculus and generalized special functions, see Kumar and Haubold 
[To] , for series expansions for all possible values of v. In many cases the nuclear factor 
S^(0) is approximately constant across the fusion window. Taking S'M(O) = for 
v = 1 and v = 2 and taking S°(0) = S(0), we obtain the extended reaction rate 
probability integral as 



7*12 



1 - ^5 12 ) 



8(g-l) 
fikT 



7T 



5(0) G\ 



3.1 

,3 



(q-l)E G ,2--L- 1 



AkT 



'0,1,1 



(2.10) 



The series representation for (]2. 10[) can be obtained as 

-S(0) 



ri 2 



1 - 7^12 ) //J.//2 



2vrr 



+ 



1 



q-1 
2^(q - 1)E G 
AkT 



8(g-l) 



7T 



fikT 
(9 - l)Bo 



1) 



r 



i 



r=0 



4A;T 

(g - i)E G 
AkT 



1F2 



q-1 

1 3 1 (q-l)E G 



q-1 2' 2' 2' AkT 



A r - In 



(g - l)E G 
AkT 



(2.11) 



where 



and 



1 
2 



1 



9-1 

-lVr 



-E>r 



1 

9-1 



+ r) + tf(2 + r) 



(§) r!(l + r)! 



(2.12) 



(2.13) 



see the Appendix for detailed evaluation. As q — y 1 in (2.10), then by using Stirlings 
formula for gamma functions, given by 



T(z + a) « (2tt)2z 



M — > 00, a is bounded, 



(2.14) 



we get the reaction rate probability integral in the Maxwell-Boltzmann 



ri2 



5i2 n x n 2 



7T/i 



5l2 ™1™2 



/i/cT 



kT 1 S(0) 



exp 
-Eg i- 



Eg 



4A;T l0 'i' 1 



d£ (2.15) 



(2.16) 



which is obtained in a series of papers by Mathai and Haubold, see for example Mathai 
and Haubold [21] • The integral in (2.15) is dominated by the minimum value of 
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^ + (^) 2 = 9(E) ( sa y)- The minimum value of the function g(E), say E , can 



be determined as 
E 



d 

dE 



+ 



kT 

and the function 
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kf 



1 I _S 

-EZEr. 2 

2 G 



0=>E, 







i fkT\ 5 
2~ 



s(£ 



G 



AkT 



36 



(2.17) 



(2.18) 



j. 

where 6 = (J^) 3 - Now, by using the Laplace method [5j |25], we can obtain an 
approximate value for (2.15) as 



ri 2 



1 



n x n 2 



8^(0)e 2 exp(-36) 



(2.19) 



In the next section we will obtain the mass, pressure and temperature for the case of 
analytic stellar models characterized by density distribution and corresponding tem- 
perature distribution considered by Haubold and Mathai [12] . 

3 Closed forms of the integral over the stellar nu- 
clear energy generation rate 

Let us consider the density distribution g(r), considered by Haubold and Mathai 
IT21 [TBI. Mathai and Haubold EU in the form 



g(r) 



Qc 



1 - 



r y 

r) 



,5>0, 



(3.1) 



where g c is the central density of the star, r is an arbitrary distance from the center 
and R is the solar radius. This density function is capable of producing different 
density distributions by choosing the free parameter 5. Now we determine the quantities 
M(r),P(r), and T(r), the mass, pressure and temperature at r. 
By the equation of the mass conservation, 



dM(r) 
dr 



4iir 2 g(r) 



we get, 



M(r) = 4vr£ c / t 



dt 



4tt 



-Q c r 



T \ 1 



5 + 3 \RJ 



(3.2) 



(3.3) 



From (3.3), we get the central density as, 

_ 3(5 + 3) M(R) 
6c ~ Att5 J&~" 



(3.4) 
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If an element of a matter at a distance r from the center of a spherical system is in 
hydrostatic equilibrium, then setting the sum of the radial forces acting on it to zero 
we obtain, 

dP(r) = Gg(r)M(r) 
dr r 2 

where G is the gravitational constant. Assuming that the pressure at the center of the 
Sun is P c and at the surface is zero, we get 



P(r) 



Pc-G 



M(t)g{t) 



dt 



47T 



Gg 2 c R 2 



t 2 

1 /r\2 



5 + 6 



2 \RJ (5 + 2)(5 + 3) \RJ 



f \ 8+2 



r \ 28+2 

2(5 + l)(5 + 3) \RJ 



Using the boundary conditions P(R) = 0, we get, P c = ^-G^g^R 2 , where 



5 + Q 



+ 



(5 + 2)(5 + 3) 2{6 + l){6 + 3)' 
By the kinetic theory of gases, for a perfect gas, the pressure is given by 

kN A 



(3.6) 



(3.7) 



P(r) 



-g{r)T{r) 



(3- 



For the temperature of interest for stellar models, we neglect the negligible radiation 
m th 

T(r) 



pressure from the total pressure and obtain from (3.8) as 

H P(r) 



kN A g(r) 

Att Gfig c R 2 
3kN A [l -(%)*] 
3 

2(6 + l)(5 + 3) vi? 



r \ 



2 \RJ 

r \ 28+2 



2 

) + 



5 + Q 



(5 + 2)(5 + 3) \R 



<5+2 



(3.9) 



-j^Gfx^ M ^ , where £ is defined in (3.7). 



The central temperature is given by T c 

Thus we have obtained the mass, pressure, and the temperature throughout of the 
nonlinear stellar model with the density distribution defined in (3.1). Next our aim is 
to obtain analytical results for stellar luminosity and neutrino emission rates for various 
stellar models. 

4 Stellar luminosity and neutrino emission rate 

The energy conservation equation states that the net increase in the rate of energy 
flux coming out of a spherical shell from inside is the same as the energy produced with 
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in the shell [H] . If we denote L r 
r, then we have 



L(r) as the energy flux through the sphere of radius 



dL r 
dr 



47rr g(r)e(r), 



(4.1) 



where e(r) is the energy produced per second by nuclear reactions of each gram of 
stellar matter. The quantity e(r) depends on the chemical composition in each gram 
of stellar matter. Here usually L r is a constant but will be equal to L at the surface of 
the star. We assume here that the star is chemically homogeneous (that is a star where 
chemical composition throughout is a constant). Also we assume the energy generation 
rate e(r) for one particular nuclear reaction. Now if we denote r 12 (^(r), T(r)) as the 
extended non-resonant thermonuclear reaction rate for the particles 1 and 2 defined by 
(2.10), then we shall consider the energy generation rate £ 12(f) which can be written in 
terms of the extended reaction rates via 



£12 [r 



g(r) 



C*g 2 G\ 



3,1 
3 



(g -!)£<, .2-1. 



4kT 



'0.5.1 



where 



C* 



Ei 2 r 12 (g(r),T(r)) 



Q 2 G{ 



3,1 



(g-l)E G | 2 -^I 

It) i 1 



4kT 



(4.2) 



(4.3) 



in which E12 is the amount of energy given off in a single reaction. It is to be noted 
that by using the asymptotic behavior of ^13( 4^^ ) [19] and as q — > 1, C* — > C, the 
composite parameter considered by [TJ, which is defined as 



C 



E 12 Ti2 

g 2 Q 2 



exp(36) 



(4.4) 



for our universe C ~ 2 x 10 4 for proton-proton fusion under typical stellar conditions 
[lj. Then from (4.1) we have the total luminosity of the star by integration as, 



L(R) 



Airr g(r)e(r)dr. 







If we are considering only one specific reaction defined in (1.1), then we have 

■H 



L l2 (R) 



4-Kr~ g(r)£\2(r)dr, 



(4.5) 



(4.6) 



where the energy generation rate is defined in (4.2) and g(r) is a suitable density 

1 1( 

(9-1)^,2--^- 



distribution realized in the sun. Writing (4.6) in terms of tvi(q(t), T(r)) we get 



L 12 (R) 



4%r z C*g z Gf 



AkT 



lo,|,i 



dr 



R 



4irr 2 E 12 ri2(g(r),T(r))dr. 



(4.7) 
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The number density rii of a particle i, for a gas of mean density g(r), can be expressed 
as 



Xi 



(4- 



where stands for Avagadro's constant, Ai the atomic mass of particle % in atomic 
mass units and Xi is the mass fraction of particle % such that ^ Xi = 1. Substituting 
ri2(g(r),T{r)) from (2.10) and using (4.8) we have 



L 12 (R) 



R 



4irr 2 E 12 [l--5 12 \N 2 A Q 2 (r 



X\X 2 
] A X A 2 



x S(0)Glj 



(g-l)7r 2 /i ( Z x Z 2 e 2 \\i-^ 
2kT(r) V h 



9- 

lo 



8(g-l) 
jikT{r) 

dr. 



71 



-1 



1 5-1 2 



(4.9) 



If we divide the "internal luminosity" L 12 (R Q ) by the amount of energy E i2 , then we 
get the total number of particles per second N 12 liberated in the reaction given by ( 1.1 ) 

as 



N- 



L 12 {R) 



12 



E 



12 



4(1- -<J 12 ] N 



X\X 2 
X ~A\A 2 



fik 



1 



x 



R r 2 g 2 (r) 31 

r^i s 

o T r a 1,3 



4(1- -5 12 ] iV : 



(g-l)TrV / Z 1 Z 2 e 2 \\2-^ 1 
2kT(r) V ^ 



: XiX 2 

[ A^~ 2 



fik 



lo,i,i 



1 



_1__3\ 

dr 



s-(o) 



5(0) 



x 



2tt« 



r( s )r( - + s ]r(i + s )r 



1 - s 



X 



(g - l)7r 2 /i / ZiZ 2 e 



2k 



,2\ 2' 



r 2 ^(r)[T(r)]-5+Mrds. 



(4.10) 



For the density distribution defined in (3.1 ) introduced by Haubold and Mathai [T3"ll2"l"] 
and the corresponding temperature distribution (3.9) we get, 

-R 



rY(r)[T(r)]~5+Mr 



o 



+ 



4ttG> 2 
3kN2 QcR . 

5 + 6 



R 

o 

f \ 5+2 





2 * 








HQ 







(5 + 2) (6 + 3) \RJ 2(5 + !)(<? + 3) \R 



r x 25+2 



dr, 
(4.11) 
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where £ is as defined in (3.7). If we put a substitution r = xR then we get 



r 2 g 2 (r)[T(r)}-^ +s dr 



g 2 c R 3 I x 2 [l-x 5 ]f 



+ 



4ttGu 2 
£ + 6 

(5 + 2)(5 + 3)" 2(5+l)((J + 3) 



-x 5+2 - 



2<5+2 



dx. (4.12) 



Putting x s = y and simplifying we obtain, 



r 2 ^ 2 (r)[T(r)]-2 + Mr 



"47rG/i£ 
3^ 



g 2 R 3 



X I yl — J/]a s [l-u(y)] s ^dy, 



(4.13) 



where w(y) is defined as 



w(y) 



5 + 6 



2 (6 + 2)(6 + 3) V + 2(5 + l)(5 + 3) 



(4.14) 



As y — >■ 0, w(y) — >■ and as y — >■ 1, w(y) — >■ 1. If we take 



v(y) 



5 + 6 



2 (5 + 2)(5 + 3) y 2(5+1) (5 + 3)' 



(4.15) 



we have v(0) = |. The minimum value of v(y) is at y = ^7"^) an d the value is 

2 ~~ 6 (3+3) (i+2)^ • Thus the minimum value is non- negative since [f^T^p decreases 
steadily from 3 to 1 for all 5 > 0. Therefore v(y) < 0. Since £ > for all 5 > 0, 
u(y) < for all 5 > 0, £ > 0. Thus [1 - M(y)] s ~i < 0. Hence < u{y) < 1 for < y < 1 
and for 5 > 0. Thus by using the binomial expansion we obtain 



oo j\ 



[i-„ (y )ri = ^^i^[„( y )] 

z — ' m! 



(4.16) 



m=0 



where (| — s) m is the Pochhammer symbol defined for a G C by 

(a)o = 1, (a) m = a(a + 1) • • • (a + m — 1), m = 1, 2, • • • ,a ^ 
T(a + m) 

r(a) ' 

whenever T(a) exists. Taking 5 = 2 we get £ = | and 



(4.17) 



(4.18) 



11 



Then from (4.13) we obtain 



R 



rY(r)[T(r)]~5+ s dr 



4vrG^ 2 



■3 £ c 2 # 3 



y |-i[l- y ]l+-i(l-I y )-ldy 



4vrG> 2 
g c K 



!5kN A 
i 



2p3 00 ('l 



^Om 1 



x / y 



;-+m-i 



2 ^ m! 2 m 

m=0 



[l-y]t+- 1 dj/ 



By using beta integral and using (4.17) we obtain 



r 2 ^(r)[r(r)]-2 +s dr 



^ , c i? 2 



15kN A 



s ~5 ,2p3 ^ r (i _ s + m ) 



' m!r(i- S ) 



X 



i r(§ + m)r(| + s) 
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reaction ( 1.1 ) as 



Now from (|4.10|) we obtain the total number of particles per second liberated in the 
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For more details on Meijer's G-function and its properties see [TJJJ [23] . Thus we 
have obtained the the total number of particles per second liberated in the reaction 
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(1.1) in terms of the density distribution considered by Haubold and Mathai [TBI 1X41 . 



5 Comparison of the pathway energy density and 
the Maxwell-Boltzmann energy density 

In Figure 1 below it can be observed that for the nuclei to react at energy E, it 

has to borrow an energy E from the thermal environment. The probability of such an 

energy is proportional to the Maxwell-Boltzmann energy exp [— The fusion will 

take place when the nuclei penetrate the Coulomb barrier keeping them apart. The 

r -i 

probability of penetration is given by the factor exp — (^) 2 • The product of these 

two factors illustrates that fusion mostly occurs in the energy window shown in the 
following Figure 1. 



1 1 



: oc exp 



■I 
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kT 



a 
a 
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□ 
□ 
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E 



8kT 10 
Energy 



Figure 1: Schematic plot of the energy-dependent factors for the reaction rate probabil- 
ity integral: Maxwell-Boltzmann energy density, non-resonant nuclear cross-section and the 
product of the Maxwell-Boltzmann density and the non-resonant cross-section. 



In Figure 2 below the pathway energy density is plotted for q = 0.7, 0.9, 1, 1.2, 1.4, 
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respectively. For different values of q we get different energy densities [Curves (a), 
(b), (c), (d), (e)]. The non- resonant cross-section is also plotted. The product of 
the pathway energy density and the non-resonant cross-section for different values of 
q namely q = 0.7,0.9, 1, 1.2, 1.4 are also plotted. It is to be noted that as q — > 1 the 
pathway energy density coincides with the Maxwell-Boltzmann energy density and also 
the fusion window for the Maxwell-Boltzmann case in Figure 1. 




Figure 2: Schematic plot of the energy-dependent factors for the extended reaction rate prob- 
ability integral: pathway energy density, non-resonant nuclear cross-section and the product 
of the pathway density and the non-resonant cross-section. 

The curves in the figure represents pathway density [l + (q — 1)^] 9-1 for (a) q = 0.7, (b) 
q = 0.9, (c) q = 1, (d) q = 1.2, and (e) q = 1.4. The curve (f) represents the non-resonant 

. The product [l + {q - 1)^] exp - 

q = 0.9, (i) q = l,Q)q= 1.2, and (k) q = 1.4. 



cross-section exp 
resented for (g) q = 0.7, (h) 



E ) 



is also rep- 
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Figure 3(a) below shows the pathway energy density defined in (2.4) for q 



1,1.2,1.3,1.4 and Figure 3(b) shows the Maxwell-Boltzmann energy density defined 



in (1.3). As q — > 1 the pathway energy density reduces to the Maxwell-Boltzmann en- 



ergy density. The pathway energy density covers many stable and unstable situations 
as the value of q varies. If the Maxwell-Boltzmann density is the equilibrium situation, 
many other non-equilibrium situations are covered by the pathway energy density. 



Pathway Energy Denisty for k=l, T=100K 



0.003 



0.001 




Maxwell-Boltzmann Energy Density 
foTk=l,T=100K 



1000 




1000 



E Energy ^ E Energy nA 

Figure 3(a): Pathway energy density for q = 1, 1.2, 1.3, 1.4 and for k = 1,T = WOK. 
(b): Maxwell-Boltzmann energy density for k = 1,T = 100K. 



6 Concluding remarks 

In this paper we have modified the energy distribution for a non-resonant reaction 
rate probability integral. The composition parameter C considered by [TJ is extended 
to C* by the pathway energy density. Analytic density distributions, considered by 
Haubold and Mathai [T3J E] , are used to obtain the stellar luminosity and the neu- 
trino emission rates and are obtained in terms of generalized special functions such as 
Meijer's G-function. The pathway energy density considered here covers many density 
functions and hence the extended reaction rate integral covers a wider class of integrals. 
Pathway energy density helps us to obtain various fusion windows by giving different 
values to q, the pathway parameter. The graphs plotted here were obtained by using 
Maple 14 in Windows XP platform. 
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A Appendix 



Series representation 



The series representation for the right hand side of (2.10) can be obtained through 



the following procedure. Here we apply residue calculus on the G-function given in 



(2.10). Consider the G-i unction 
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The right hand side is the sum of the residues of the integrand. The poles of the gam- 
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mas in the integral representation in (2.10) £1X6 clS follows. 
Poles of T(s) : s = 0, -1, -2, . . . ; 

Poles of T (| + s) : s = -§,-§,-§,... ; 
Poles of r(l + s) : s = -1, -2, -3, . . . . 

Here the poles of T(s) and r(l + s) will coincide each other at all points except at s — 0. 
Note that the pole s = is a pole of order 1: s — — |, — |, — |, . . . are each of order 1: 
s = —1, —2, —3, . . . are each of order 2. We know that 
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1 if r = 0, 

The sum of the residues corresponding to the poles s = is given by 
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where iF% is the hypergeometric function defined by 

(a L x 



iF 2 (a; b,c; x) = ^ 



r=0 



(b) r (c) r r! 



To obtain the sum of the residues corresponding to poles s = —1, —2, —3, .... of 
order 2, we proceed as follows: 
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where ^{z) is a Psi function or digamma function (see Mathai (T9J and ^(1) 
0.5772156649 ... is Euler's constant. Now 
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Then by using (A.7), (A. 8 



) and (A. 9) we get, 
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Thus from (A.5), (A.6) and (A.10) we get (2.11). 
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